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Abstract 

Current submillimetre surveys are hindered in their 
ability to reveal detailed information on the epoch of galaxy 
formation and the evolutionary history of a high-redshift 
starburst galaxy population. The difficulties are due to the 
small primary apertures (D < 15-m) of existing submil- 
limetre telescopes and the limited sensitivities of their first 
generation of bolometer cameras. This situation is chang- 
ing rapidly due to a variety of powerful new ground-based, 
airborne and satellite FIR to millimetre wavelength facil- 
ities. Improving our understanding of the luminosity and 
clustering evolution provides the motivation for conduct- 
ing cosmological submillimetre and millimetre surveys. It 
is therefore important that we quantify the limitations of 
the future surveys and the significance of the results that 
can be drawn from them. In this paper we present sim- 
ulated surveys which are made as realistic as possible in 
order to address some key issues confronting existing and 
forthcoming surveys. We discuss the results from simu- 
lations with a range of wavelengths (200/im - 1.1mm), 
spatial resolutions (6"- 27") and flux densities (0.01-310 
mJy). We address how the measured source-counts could 
be affected by resolution and confusion, by the survey sen- 
sitivity and noise, and by the sampling variance due to 
clustering and shot-noise. 
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1. Introduction 

By early 2001 the initial extensive programme of extra- 
galactic SCUBA (850/im) surveys conducted on the 15- 
m JCMT will be completed. The targets of these sub- 
mm surveys include lensing clusters (Smail et al. 1997), 
the Hubble Deep Field (Hughes et al. 1998), the Hawaii 
Deep Survey fields (Barger et al. 1998, 1999), the Canada- 
France Redshift Survey fields (Eales et al. 1999, Lilly et al. 
1999) and the UK 8 mJy SCUBA surveys of the ISOPHOT 
ELAIS and Lockman Hole field (in prep.) These blank- 
field sub- mm surveys, which cover areas of ~ 0.002 — 
0.2 deg 2 , allow a measure of the 850/im source-counts be- 
tween flux densities of l-15mJy. The relatively small sur- 
vey areas are due to the low mapping speed (a combina- 



tion of the field-of-view and sensitivity of the bolometer ar- 
ray, and the necessary overheads to fully sample the focal- 
plane with > 0.5FA feedhorns, e.g. 16 separate secondary- 
mirror positions are required to fully sample the SCUBA 
array). The restricted 850/im flux range in the measured 
source-counts is the result of a source-confusion limit of 
~ 2 mJy (due to a small primary aperture or map reso- 
lution) at faint flux densities and the small survey area 
at the bright end, given the intrinsic low source-density 
of bright sub-mm sources (e.g. iV[585o Mm > 20mJy] < 
70 sources/deg 2 ). 

The limitations on the accuracy of our understanding 
of high-z galaxy evolution from sub-mm surveys has been 
described elsewhere (Hughes 2000). To improve the con- 
straints on the competing evolutionary models, future sub- 
mm and mm surveys must extend their wavelength cover- 
age, increase their survey area and sensitivity. To quantify 
the advantages of these forthcoming surveys we have sim- 
ulated the extra-galactic sky at various sub-mm and mm 
wavelengths (200-3000/im), spatial resolutions (Ofwhm ~ 
1 — 120") and sensitivities (S > 0.01 mJy). We include in 
the simulations contributions from CMB primary fluctu- 
ations, S-Z clusters, an evolving extra-galactic population 
of star-burst galaxies and high-latitude galactic-plane cir- 
rus emission. These simulated surveys allow us to compare 
the measured source-counts for a wide variety of current 
and future FIR and sub- mm/mm telescopes (e.g. SIRTF, 
FIRST, MAP, PLANCK, BLAST, SOFIA, CSO, JCMT, 
LMT, GBT, ALMA). 

We discuss these multi- wavelength simulations in more 
detail elsewhere (Hughes & Gaztahaga 2000). In this pa- 
per we concentrate on illustrating the effects of galaxy 
clustering and low spatial resolution on the measured source- 
counts at 1100/im, 850/im, 350/im and 200/im. 

2. Simulations 

The starting point for our simulations is a P 3 M N-body 
simulation with 200 3 particles in a 600/i _1 Mpc box that 
has the same matter power spectrum as that measured for 
APM Survey galaxies (see Gaztahaga & Baugh 1998, and 
references therein). Thus, the 2-point correlation in the 
(dark) matter particles is identical to the one measured in 
the local (z ~ 0.15) galaxy universe. Using a single N-body 
output (as opposed to a light-cone output) corresponds to 
having the matter clustering pattern fixed in comoving 
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coordinates (stable clustering). If redshift evolution of the 
two-point correlation function £2(^2) is parametrized as: 

Ur,z) = {l + z)-^ Ur), (1) 

then, for stable clustering we have e ~ at small scales 
{e.g. see Gaztahaga 1995). This has less clustering evo- 
lution than pure matter gravitational clustering {e.g. lin- 
ear or non-linear growth, where e ~ 1) and might de- 
scribe some models in which galaxies are identified with 
high-density matter peaks (peaks move less than particles, 
which results in less evolution). This might be adequate 
for star-forming galaxies if we are selecting rarer (more 
massive) objects as we go to higher redshifts. It is in agree- 
ment with the strong clustering observed in Lyman-break 
galaxies at z ~ 3, which is comparable to the clustering 
of present-day galaxies (see Giavalisco et al. 1998). 

We next implement a model for the evolution of the 
galaxy population and set the geometry of a light-cone 
observation. To relate spatial coordinates r and observed 
redshifts z in the simulation output we use a spatially 
flat metric with the critical density Q = 1 and zero cos- 
mological constant (other cosmological parameters will be 
considered elsewhere). We transform the N-body matter 
simulation into a mock galaxy catalogue of angular and 
redshift positions by the following steps: 

— i) select an arbitrary point in the simulated box to be 
the local 'observer'; 

— ii) apply a given survey angular mask {e.g. a square 
with 1 degree on a side); 

— iii) include in the mock catalogue a simulated particle 
at redshift z from the observer with probability given 
by some selection function ip{z); 

— iv) assign a luminosity L = L{\) to this particle ac- 
cording to the luminosity function 4>[L,z] for each one 
of the observer filters A, at the corresponding particle 
redshift z. 

We replicate the N-body (periodic) simulation box to 
cover the total extent of the survey (up to z = 6, beyond 
which we impose an exponential cut-off in the expected 
number of galaxies). By comparing the results from dif- 
ferent box sizes we have verified that this replication of 
the box does not introduce any spurious correlations on 
large scales. The selection function ijj{z) is the normalized 
probability that a galaxy at redshift z is included in the 
catalogue, and is proportional to the estimated number of 
galaxies at this coordinate: 

rL 2 (z) 

ip{z) oc / dL (j>[L,z\. (2) 

where 4>[L,z] is the luminosity function and L\{z) and 
L<i{z) are the luminosities corresponding to the lower and 
upper limits in the range of apparent fluxes used to build 
the galaxy sample or catalog under study. Thus our pre- 
scription for galaxy formation is that, on average, the 
probability to find a galaxy is simply proportional to the 
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Observed Wavelength (yum) 

Figure 1. Radio to IR SEDs of low-z starbursts, quasars, 
ULIRGS, Seyferts normalised at 60fim. The curves repre- 
sent fits to the SEDs of Arp220 (solid), M82 (dashed) and 
Mkn231 (dashed- dotted). Since the rest-frame SEDs show 
a large dispersion in the spectral slope at A < 60 jam then 
only the number-counts of the highest-z sources (z > 3) in 
the shortest wavelength simulations (200-300jim) are sen- 
sitive to the exact choice of galaxy SED type or its redshift 
evolution. 

(dark) matter densities, so that the number density of 
galaxies automatically determines how massive are the 
matter peaks that our galaxies are tracing. Other "bi- 
asing" schemes can also be implemented, but note that 
different models for galaxy formation can give different 
clustering evolutions for a given (f)[L,z], so that both of 
these evolutionary patterns need to be better constrained 
by the future data if we want more realistic simulations. 

To include the star-formation history of the galaxies 
in our simulations we have selected models describing the 
SEDs shown in Figure [l| to extrapolate the local IRAS 
60/im local luminosity function <j)[L,z = 0] = cj)[L/L*] 
(Saunders et al. 1990) to longer wavelengths. We have in- 
troduced an exponential cutoff for objects fainter than 
L = 1O 8 L and adopted a model of pure luminosity evo- 
lution, (j)[L,z\ = cj)[L/L*], with the following redshift de- 
pendence: L*(z) = (1 + zf L*{z = 0) for < z < 2.2; 
L*{z) = L*{z = 2.2) for 2.2 < z < 6; an exponential cutoff 
at z > 6. 

In this paper we use a single SED model of Arp220 for 
all galaxies in our sample. As will be shown below (Figure 
[|), using a different SED makes no significant difference 
in the number counts or redshift distributions even for 
bright sources at the shortest submillimetre wavelengths. 
We have also assumed no evolution for the SED: if the 
evolution occurs within the spread in the models shown 
in Figure |l|, this has little effect for A - 200 - 1100/xm. 

Thus, after these transformations we end up with a 
galaxy catalogue of angular positions, redshifts and lumi- 
nosities, for each of the observer filters A. We next smooth 
the resulting projected distribution with a spatial resolu- 
tion appropriate to the observational configuration we are 
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Figure 2. A confusion-limited (la = 0.5 mJy) JCMT 1 deg 2 area simulated survey at 850jim with Ofwhm = 15". The 
large and small boxes show the 0.1 deg 2 area and the HDF-like regions shown in Figures & respectively. 



simulating: Ofwhm — 1-2 A/D where D is the telescope 
aperture. Finally we add white noise M to the angular 
catalogue to emulate the signal-to-noise S/M appropriate 
for the integration time, t, the instrumental sensitivity and 
observational conditions (Noise Equivalent Flux Density, 
NEFD) of the mock survey we want to emulate. 

We describe here the results from simulated surveys 
at 1.1 mm, 850/im, 350/im and 200/im with spatial resolu- 
tions of 6", 15", 9"and 25 // corresponding to surveys on the 



50-m LMT ( |http://binizaa.inaoep.m^), 15-m JCMT, 10- 
m CSO and 2-m BLAST ( http: / / www.hep.upenn.edu/blast / ) 
respectively. The full simulations have a flux dynamic range 
of S\ = 0.01 — 310 mJy at all wavelengths and cover an 
area of 1 sq. degree. Subsets have been extracted to deter- 
mine the source-counts from surveys comparable in area 
(6-400 arcmin 2 , Table g - survey A and B) to the SCUBA 
surveys of the Hubble Deep Field and Hawaii Deep Fields, 
the lensing cluster survey, Canada-France Redshift survey 
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LMT Surveys 


NEFD = 


10 mJy sec 1/2 , n ~ 100/rni 


NEFD 4 mJy sec 1/2 , n - 70/im 


Survey 


Area (sq. arcmins) 


3-cr depth 


N (5 — a galaxies) 


3-cr depth N (5-cr galaxies) 


A 


6 


0.1 mJy 


18 


0.05 mJy 33 


B 


350 


0.9 mJy 


204 


0.4 mJy 515 


C 


1000 


1.6 mJy 


320 


0.6 mJy 1100 


D 


3600 


3.0 mJy 


400 


1.2 mJy 1600 



Table 1. Number of 5-a galaxies detected in alternative 50-hour BOLOCAM 1.1 mm surveys with a conservative Noise 
Equivalent Flux Density (NEFD) of NEFD = 10 mJy sec 1 ! 2 a i 11 mm, during the initial commissioning phase (e.g. 
with a primary aperture r.m.s. surface accuracy of n ~ lOO/am r.m.s.) of the LMT; and later during routine operation 
with an improved surface accuracy (of r) ~ 70/^m) and 1.1mm NEFD = 4 mJy sec 1 ! 2 . BOLOCAM sensitivities allow 
for overheads. Surveys A and B represent LMT surveys similar in area to the SCUBA survey of the HDF and the UK 
8 mJy wide- area survey respectively. Surveys C and D represent future wide- area LMT surveys. 
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FigureS. Simulated surveys covering 0.1 deg 2 : A confusion-limited (3a = 0.05 mJy) LMT survey at HOOfim with 
Ofwhm = 6" (top-left); a Ofwhm = 21" low resolution HOOfim survey (top-right); a Ofwhm = 15" deep (la = 
0.5 mJy) 850fim survey (middle-left); a Ofwhm = 15" shallow (la = 2.5 mJy) 850fim survey (middle-right); a 
confusion-limited (Ofwhm = 9") 350jim surve y (bottom-left); a confusion-limited (Ofwhm = 25") 200fim survey 
(bottom-right). These figures look much better in http:/ /www. inaoep.mx^gazta/ simsubmm.pt 
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fields and the UK 8mJy ELAIS and Lockman Hole sur- 
vey. Additionally larger-area surveys (0.3-1.0 deg 2 , Table 
1 - surveys C and D) more appropriate to the LMT and 
BLAST are also considered. 

Figure|| shows a large projected (1 deg 2 ) map at 850/im 
with some identified sub-regions within which we perform 
the source-counts analysis. Although there is significant 
clustering in this figure, the overall distribution is already 
quite homogeneous. Figure || shows a comparison of maps 
at different resolutions and wavelengths from the 0.1 deg 2 
region in the center of the larger map. 

In general, one can recognise the brightest sources in 
an individual map at any other wavelength. Consequently 
it is possible to derive a robust constraint on the redshift 
of an individual galaxy from the relative intensities at dif- 
ferent wavelengths (e.g. Hughes 2000). Nevertheless, even 
if we choose the sensitivity limits at each wavelength to 
have similar mean redshift s, the redshift distribution for 
the 200/im galaxies is quite different from that of 850/im 
(e.g. see Figure |]). Therefore care must be taken in de- 
termining the appropriate depths of the complementary 
surveys when the goal is to measure millimetre to submil- 
limetre colours to estimate the redshift distribution. These 
redshift issues will be quantified in more detail elsewhere 
(Hughes & Gaztahaga 2000). 

3. Source-Count Analysis 

To demonstrate the simulations accurately reproduce the 
input models we show in Figure || the match between the 
expected number of galaxies at 850/im and 200/im in a 
given log redshift bin, N(z), as given by the input se- 
lection function, and the measured counts in the redshift 
simulations (before including the noise or the finite res- 
olution to the map). Figure |] also shows a comparison 
of the results for the brighter counts using two different 
SEDs models: Arp220 and M82. As shown in Figure ||, 
the difference in the SED models is larger at A > 60/im, 
so that for mean redshifts of z c± 2 — 3, the choice and/or 
evolution of the galaxy SED will only affect the counts at 
the shortest sub-mm wavelengths (e.g. 300 — 200/im). The 
right-hand panel of Figure ^ illustrates that at 200/im the 
effects of using different SED models on the bright counts 
(~ 100 mJy) are small. The results at fainter 200/im fluxes 
or at other wavelengths are not shown as the differences 
in the histograms are indistinguishable for different SEDs. 
However the effect on the counts of the brightest 200/im 
sources is also negligible compared to the sampling effects 



due to clustering discussed in § p.3 



The number counts and photometry of the sources in 
our simulated angular maps were determined using SEX- 
TRACTOR (within the STARLINK package GAIA). In 
Figure || we compare the predicted and measured number- 
counts of the sources extracted from the full 1-deg 2 sim- 
ulations at 1100, 850 and 200/im. Again there is excellent 
agreement between the simulated source-counts and the 



model down to the confusion flux limit of the different 
surveys at their respective resolutions. This illustrates the 
obvious point that provided the survey is of sufficient area 
and sensitivity then neither resolution, projection or noise 
are very important in extracting counts for objects above 
the confusion-limit. 

3.1. Resolution and confusion 

Below the confusion limit the counts flatten as faint sources 
merge to form brighter objects (see Figure ||). For ex- 
ample, at 850/im the counts are affected by confusion at 
'S'sso^m < 2mJy, whilst at 1.1mm with G^resolution (e.g. 
LMT) the counts are still unaffected at Si.i mm ~ 0.1 mJy. 
Even at lower resolution (27") the measured 1.1 mm source- 
counts recover the input model down to 5i.i mm ~ 3mJy 
due to the low density of intrinsically luminous sources (> 
lO 12 I/0). Similarly the 200/im counts recover the model 
down to a confusion limit of 18 mJy at 25 // resolution. Hence, 
even at resolutions of ~ 30" the overlapping of individual 
source PSFs is negligible, regardless of whether it's due 
to the random line-of-sight projection of galaxies at dif- 
ferent redshifts, or because of a high-amplitude of clus- 
tering. The significance of this result is that it is there- 
fore possible to combine large-aperture (50-m LMT, 15- 
m JCMT) 3000-450/im wavelength surveys and relatively 
small primary aperture (e.g. ~ 2.0m BLAST) sub-mm 
and FIR surveys (500 — 200/im) with very different resolu- 
tions (0FiR/0 mm ~ 5) and still derive meaningful colours 
of sub-mm-selected galaxies, and hence constrain their 
redshifts (Hughes 2000). 

3.1.1. Extra-galactic Background 

The 1.1mm extra-galactic background with flux S > 0.5 
mJy due to the discrete sources in our simulations is 1.31 x 
10 -10 Wm -2 sr -1 , in excellent agreement with FIRAS (Far 
Infrared Absolute Spectrophotometer) results from COBE 
residual measurement of the extra-galactic far infrared 
background (e.g. Dwek et al. 1998 and references therein). 
This implies that surveys conducted at the full spatial res- 
olution of the LMT could resolve ~ 100% of the mm back- 
ground, whilst the deepest SCUBA surveys todate have 
resolved only 30-50% of the sub-mm background. 

3.2. Depth and noise 

In the previous subsection we indicated that provided sur- 
veys have sufficient sensitivity that they can reach their 
confusion limit then accurate source-counts can be mea- 
sured from large- area surveys. However the simulations 
also illustrate the way in which the counts are severely 
affected for shallow wide-area surveys with a flux-limit 
well above the confusion limit when the number-counts are 
steep. For example Figure || shows the simulated counts 
from a 0.2 deg 2 850/im SCUBA survey with a 1-a sensitiv- 
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redshift z redshift z 

Figure^- Redshift distribution, N{z), of simulated sources (histograms) compared to input model (solid-lines). Left- 
panel shows N(z) at 850fim for sources brighter than 0.01, 1.0, 10.0 mJy. Right-panel compares N(z) at 200fim for 
sources brighter than 10 and 100 mJy: The dotted and continuous lines correspond to the use of Arp220 and M82 as 
the template to model the galaxy SED. 
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Figure 5. A comparison of the input model (solid-lines) 
and the extracted source- counts from the simulated 1-deg 2 
surveys at 200fim, 850fim and 1.1 mm. Note how the sim- 
ulated counts are depleted with respect to the input model 
at fluxes fainter than the confusion limit, which is set by 
the spatial resolution of the survey. The 1.1 mm counts at 
different resolutions (6" and 21" ) are also shown. 



ity of 2.5 mJy. The extracted counts accurately reproduce 
the model down to a flux density of ~ 12 mJy, below which 
a steep increase in the counts is observed due entirely to 
the noisy background. Hence individual sources identified 
in shallow sub-mm surveys at the 3-5a level may be spu- 
rious. These additional sources can have a non- negligible 
positive contribution to the overall number counts (com- 
pared to the situation in §3.1 where the counts flatten at 
the confusion-limit of the survey). This is due to the steep 
decrease in the galaxy counts which results in an increas- 
ing contribution of the noise counts as the surveys move 
to brighter flux limits. As the noise spectrum is well de- 
termined (e.g. dashed-line in Figure ||) a correction can be 
applied to remove the excess counts (Hughes & Gaztahaga 
2000). 

3.3. Clustering and shot-noise 

Another important aspect illustrated by the simulations is 
the effect induced on the counts by the sampling variance 
of the large-scale galaxy clustering. This is most signifi- 
cant in the deepest surveys which necessarily provide the 
smallest survey areas (< 0.1 deg 2 ). In Figure ff] we compare 
3 surveys identical in area (~ 6arcmin 2 ) to the SCUBA 
surveys of the Hubble Deep Fields (Hughes et al. 1998) 
and similar to those of the Hawaii Deep Fields (Barger et 
al. 1998) and lensing galaxy clusters (Smail et al. 1997). 
We find a factor of 3-10 variation in the extracted counts 
from these deep confusion-limited surveys. The three sim- 
ulated deep surveys and the original SCUBA HDF survey 
are shown in Figure |[ For the SCUBA surveys we would 
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Figure 6. Extracted number-counts at 850/nm from two 
(3a ~ SmJy) 0.2 deg 2 simulated surveys (dashed-lines). 
The symbols show the actual observed SCUBA 850\im 
number-counts from the surveys described in §1 and the 
solid line correspond to the the model described in §2 
(which is the input to the simulations). The counts from 
the noise map (la = 2.5 mJy), which was added to the 
raw simulation, are shown as a dotted line. Because these 
are not deep surveys, the noise dominates the counts pro- 
ducing spurious simulated sources at 8-12 mJy. This has 
the opposite effect to the confusion produced by the finite 
resolution shown in Figure 0. 



need an area over 100 times larger if we want to reduce 
the variance in the counts to the few percent level. Note 
that besides the intrinsic clustering, shot-noise plays an 
important role in this variance, at least for the brighter 
(less numerous) sources. A detailed analysis of this will be 
presented elsewhere (Hughes & Gaztahaga 2000). 



4. Conclusions and Future Work 

We have presented simulated surveys which are made as 
realistic as possible in order to address some key issues 
confronting existing and forthcoming millimetre and sub- 
millimetre surveys. These simulations assume a model for 
the luminosity and clustering evolution which is poorly 
constrained by current data. Improving our understand- 
ing of this evolution provides the motivation for conduct- 
ing submillimetre and millimetre surveys. It is therefore 
important that we quantify the limitations of these sur- 
veys and the significance of the results drawn from them. 




flux density (mJy) 

Figure 7. Histograms representing the 850jim source- 
counts from three different small area surveys (~ 
6 arcmin 2 ) as shown by the small boxes in Figure ||. These 
are similar in area and depth to SCUBA 850\im survey of 
the HDF (see Figure The solid-line and data are the 
same as in Figure |^. 




Figure 8. Three deep 850 fim surveys extracted from the 
simulation in Figure || with same area and depth as the 
SCUBA 850jim survey of the HDF (shown as a circle in- 
side a black square for comparison). 
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In this paper we discuss the results from survey sim- 
ulations with a range of wavelengths (200/im - 1.1mm), 
spatial resolutions (6"- 27") and flux densities (0.01-310 
mJy) . Several issues have been addressed related to the es- 
timation of the measured source-counts in future surveys: 
resolution and confusion; survey sensitivity and noise; and 
sampling variance due to clustering and shot-noise. 

The redshift information, determined from the mil- 
limetre and submillimetre colours of individual sources 
(Hughes 2000), is an essential ingredient to discriminate 
between the possible evolutionary models that describe 
the star formation history of galaxies. 

Clustering in future wide-area submillimetre surveys 
might also help to understand some of these issues. Both 
the number counts and the two-point clustering pattern 
depend strongly on the cosmological model, but its combi- 
nation can break some degeneracies. For example if galaxy 
formation occurs in the rare peaks of the underlying mat- 
ter field then the resulting galaxy clustering should be 
stronger than if galaxy formation occurred in random places, 
The measure of the coherence of the clustering pattern, 
as described by the higher-order statistics (which is ap- 
parent in our simulated maps), can also be used to test 
how star forming galaxies trace the underlying mass (e.g. 
Gaztahaga 2000, and references therein). This issues will 
be address in more detail elsewhere. 



References 

Barger, A.J. et al. 1998, Nature, 394, 248 
Barger, A.J., Cowie, L.L., Sanders, D.B. 1999, ApJ, 518, L5 
Dwek, E. et al. 1998, ApJ 508, 106 
Eales, S.A. et al. 1999, ApJ, 515, 518 
Gaztanaga, E. & Baugh, C. 1998, MNRAS, 294, 229 
Gaztanaga, E., 1995, ApJ, 454, 561 



Gaztanaga, E., 2000, |astro-h/0003160| 



Giavalisco, M. et al. 1998, ApJ, 503, 543 
Hughes, D.H. & Gaztanaga, E., 2000, MNRAS, in preparation 
Hughes, D.H. et al. 1998, Nature, 394 241 



Hughes, D.H., 2000, astro-ph/0003414 



Lilly, S. et al. 1999, ApJ, 518, 441 
Saunders, W. et al. 1990, MNRAS, 242, 318 
Smail, I. et al. 1997, MNRAS, 490, L5 



